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Abstract 

We study theoretically, the photoluminescence properties of a single quantum dot in a microcav- 
ity under incoherent excitation. We propose a microscopic quantum statistical approach providing 
a Lindblad (thus completely positive) description of pumping and decay mechanisms of the quan- 
tum dot and of the cavity mode. Our analytical results show that strong coupling (SC) and 
linewidths are largely independent on the pumping intensity (until saturation effects come into 
play), in contrast to previous theoretical findings. We shall show the reliable predicting character 
of our theoretical framework in the analysis of various recent experiments. 

PACS numbers: 
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Cavity quantum electrodynamics (CQED) studies the interaction between a quantum 
emitter and a single radiation-field mode. When an atom is strongly coupled to a cavity 
mode it is possible to realize important quantum information processing tasks, such as 
controlled coherent coupling and entanglement of distinguishable quantum systems 2|. In 
this respect solid-state devices, and in particular semiconductor quantum dots (QD), are 
the most promising architectures for the possibility of miniaturization, electrical injection, 
control and scalability. Indeed semiconductor quantum dots provide nanoscale electronic 
confinement resulting in discrete energy levels, and an atom-like light-matter interaction. 
Cavity quantum electrodynamics addresses properties of atomlike emitters in cavities and 
can be divided into a weak and a strong coupling regimes. For weak coupling, the spon- 
taneous emission can be enhanced or reduced compared with its vacuum level by tuning 
discrete cavity modes in and out of resonance with the emitter. However, when the interac- 
tion strength overcomes losses, the system enters the so-called strong coupling (SC) regime 
[sl, Q]. In this case the usual irreversible spontaneous emission dynamics changes into a 
reversible exchange of energy between the emitter and the cavity mode. 

Commonly, these solid state systems in the SC regime are characterized with respect to 
their behaviour under optical incoherent pump excitation . In this situation, 

the pump creates an incoherent population of electron-hole pairs typically far above reso- 
nance which relaxes incoherently into the QD through various different mechanisms before 
being emitted by recombination. Moreover, some experimental observations of SC of QDs 
in microcavities, indicate that the cavity mode is weakly coupled to others various electronic 
transitions of the system which would contribute in feeding the cavity mode as well, thus 
resulting in a second incoherent pumping channel. The specific interplay of photon and 
exciton pumping and decay results in a mixed quantum steady state that infiuences consid- 
erably the observed spectra as we are going to show in detail (see Fig. lb). An analogous 
incoherent excitation would also be achieved in the case of current injection, highly desired 
for the development of optoelectronic quantum devices. Hence, appropriate theoretical mod- 
eling of SC with semiconductor QDs under incoherent excitation with a reliable predicting 
character are sought in order to fulfill the great expectations nowadays attended from future 
implementations of SC in QD systems. Of particular interest is the analysis for increasing 
pump excitation while the system is brought into the nonlinear regime. 

Recently, Laussy et al. [^ reports on the first theoretical analysis of this situation. Their 
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master equation includes pumping as absorption from incoherently populated reservoirs and 
spontaneous emission of excitons (radiative recombination) as well as of photons (cavity 
losses). Their bosonic model is well suited for the investigation of the very low excitation 
regime where linear optics dominates. It is currently heavily exploited (see e.g. Refs. [s], 
^) to analyze the SC of single QDs as it, at least at first sight, displays an impressive 
predicting character maintaining a relative easy picture and analytical results. Nevertheless 
the results obtained within their approach display some puzzling features which deserve 
careful investigations. As instance, the dependence of the polariton broadenings, of the 
Rabi splitting and of the spectra on the incoherent pumping rates as well as amplification 
effects are commonly unexpected results of an harmonic linear dynamics. 

The authors have the merit to have pointed out clearly the importance of the effective 
feeding the other electronic transitions (viz. the other QDs existing in the sample), weakly 
coupled to the cavity mode. This scenario sets a definite departure from atom-QED models 
with expected new peculiar features (see e.g. Fig. [lb and Fig. [lb) which are expected to 
influence even the quantum statistics of the emitted quanta. Anyway, although having a 
clear and reasonable physical meaning, the master equation of Ref. [7| is not of Lindblad 
form and then there is no guarantee that it safely generates a completely positive (CP) open 
dynamics Q]. Actually the obtained photon and exciton populations (see e.g. Eq. (3) of Ref. 
7|) can diverge as soon as pumping rates approaches the values of the dampings and even 
assume negative values as well as it starts to overcome. This highly undesired feature puts 
into question also the results obtained at lower pumping rates when populations remain 
finite and positive. The puzzling features, mentioned above, could be the tail (at lower 
excitation densities) of such problems. Moreover, the need for a fresh, flexible and direct 
tool for the analysis of state-of-the-art experiments is leading various groups to take 

advantage of this machinery. These possible artifacts could affect strongly the interpretations 
of the experimental data and then have negative repercussions on applications and device 
modeling. 

In this letter we provide a theoretical model able to describe SC of a single QD in a semi- 
conductor microcavity under incoherent excitation in the low and intermediate excitation 
regimes where nonlinear optical effects start to appear. In order to avoid inconsistencies 
and artifacts, we will start from a microscopic (though simple) description of the pumping 
and decay mechanisms of the QD and the cavity mode. The picture that we have in mind is 
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that of two strongly coupled subsystems (the cavity mode and the quantum dot excitations), 
each in interaction with two independent reservoirs providing both damping and pumping 
mechanisms. The laser generates electron-hole pairs in the continuum wetting layer which, 
subsequently, relax into the dot (the pump contribution) by means of incoherent scat- 
tering mechanisms such as exciton-exciton scattering mediated by phonons. At the same 
time, the cavity mode is weakly coupled to others various electronic transitions of the system 
which would contribute in feeding the cavity mode as well, resulting in a second incoherent 
pumping channel (Pa)- In the limit of weak excitation density the resulting dynamics co- 
incides with that of two strongly coupled harmonic oscillators in the presence of reservoirs. 
Nevertheless the obtained analytical results significantly differs from those of Ref. [3]. In 
particular, in the low excitation regime we obtain a very simple analytical expression for the 
PL spectrum which is definite positive, do not display any amplification effect as well as any 
change of the Rabi splitting and of the broadenings as a function of the pump intensities. 
For higher pump densities saturation effects and even lasing effects come into play. 
The master equation for this strongly interacting system can be written as 

p = t[p,Hs]+C^MC + ^QD, (1) 

where the system Hamiltonian reads 

Hs = uJaO^a + cr+a^ + g{a) + a a^) , (2) 

with g being the interaction strength between the cavity mode (with annihilation operator a) 
at energy uja and the lowest energy {u^) quantum dot exciton with transition operator from 
the ground state to the exciton level 0"+ = \e) {g\ ((t_ = (t^). The superoperators C^,ic ^^"^ 
Cqj^ describe the interaction of the cavity mode and of the QD with the reservoirs providing 
both damping and pumping mechanisms. Transmission and diffraction losses of the cavity 
mode can be modeled, within a quasimode picture, as an effective coupling (7^^) with an 
ensemble of electromagnetic modes through the output mirror [9|, [lOj . When dealing with 
optical frequencies it can safely be regarded as a zero temperature reservoir [l^ . The second 
mechanism describes the incoherent optical pumping of the cavity. It takes into account that 
in these samples there may be QDs or more generally electronic transitions weakly coupled to 
the cavity, in addition to the one that undergoes SC, determining an effective pumping of the 
cavity mode {Pa)- By applying the usual Born-Markov and rotating- wave approximations. 



we obtain, 

^Mc — " ^ ^" {2apa) — a)ap — pa) a) + -^{2a) pa — aa^ p — paa^) , (3) 

where 7a = 7^^ + 7^ contains contributions from both the reservoirs and Pa is the total 
pumping rate depending on the populations of the electronic levels weakly coupled to the 
cavity mode. Assuming only direct and weakly pumped electronic transitions we obtain 
Pa = 7^ = where (rii) is the population of the i-th level and = 

'^liOililf + il^i ~ ^aY], being 7^ the inverse of the dephasing time of the z-th transition. 
The form of the two terms in eq. ([3]) clearly shows the physical emerging picture different 
from that adressed in Ref. [3]. Indeed, the cavity experiences radiative spontaneous emission 
guided by vacuum fluctuations (proportional to 7^^), while the second reservoir provides both 
emission (proportional to 7f ((?^i) + 1)) and absorption (to 7f (nj)). 

The material excitation strongly coupled to the cavity mode is also under the influence of 
two different reservoirs: Cqj^ = Cq^ + Cq^. The first term describes spontaneous emission 
in all the available light modes except the cavity one, 

^QD = (7x/2)(2(T_p(T+ - (T+(T_p - pa+a_) . (4) 

The latter 

^QD = (ScTeiPCTje - ajjp - pajj) , (5) 

j 

is responsible for the incoherent excitation of the QD (e.g. via phonon induced relaxations) 
from the j-th levels at higher energy that get populated by optical pumping. The aj are 
temperature-dependent phonon-assisted scattering rates from the j-th excitonic levels to the 



state |e) The transition from the lowest exciton state |e) towards higher energy states 
can be safely neglected since the experiments of interests are carried out at KT << Uj — uj^- 
This master equation induces an open hierarchy of dynamical equations. In the SC regime 
the only meaningful truncation scheme is the one based on the smallness of the excitation 
density which allows to truncate with respect to the number of photon number-states to be 
included. The inclusion of one-photon states only gives linear optical Bose-like dynamics. 
On the other hand, once also 2-photon states are taken into account, saturation and lowest 
order (x*-^"*) nonlinear optical effects arise. Within our model we are able to go beyond 2- 
photon nonlinear dynamics. We shall show essentially non-perturbative calculations, in the 
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sense that with respect to the chose pumping intensity, we shall always check the convergence 
of the presented results to higher-order number-state truncation, see e.g. Fig. [1^. 

In the weak excitation regime (linear dynamics) {{(Jee) << {(^gg) — 1); from dt (A) = 
Ti {Ap), we obtain 

dt {a) = -iua (a) - ig {a_) 
dt (o-_) = -ioj^ (cr_) - ig (a) , (6) 

which is a system of two coupled damped oscillators in the absence of any external driving, 
u)c = ct^c — «7c/2 with c = a,x . We also observe that the dephasing rates do not depend on 
the pumping populations (rij) in contrast to the results of Ref. [3]. We are interested in cal- 
culating the steady-state emission spectrum: S{uj) = limt_,oo 2Re (^a^(t)a{t + r)) e^'^'^dr. 
According to the quantum regression theorem, two-time correlations {An{t)Am{t + t)) fol- 
low the same dynamics of one-body correlation functions {Amir)) but with the one-time 
correlation {An{t)Am{t)) as initial conditions. In our specific case the initial conditions are 
provided by the steady-state cavity occupation = limj^oo (c^^cf) and by C = limt^oo{0'^<^+) , 
We obtain. 

Pa , 9^ {la + lx){laPx + IxPa) /^x 

7a 7a g^ila + IxY + 7a7x \^a " ^x\ 

C = , ^ ^ {ria - n^) , (8) 

^a — ^x 

where rix can be obtained from Eq. (I7j) just exchanging the labels a and x, and = 
We end up with the following very simple expression, for the PL spectrum 

Sico) = 2Re \^tz^IhL±lB , (9) 

where the complex polariton energies determining the spectrum resonances are given by 

Rii,2) = ^^^^^T\^Ag' + {Co^~Co^Y. (10) 

The analytical structure of eq. ([9]) shows peculiar non-standard features typical of this 
system, see Fig. [T] Indeed, in addition to the usual dependence in the denominator on 
the two Rabi frequencies - responsible for the presence of the double-peak structure - we 
can appreciate a numerator carrying a resonant term proportional to the difference between 
photon and exciton occupations plus another term proportional to the sole photon steady- 
state density Ua- Both terms modulate the Rabi resonances with respect to the pumping 
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FIG. 1: (Color online), (a): PL spectrum (continuous line) calculated from Eq. ([Oj). We used 
A = LVa - iVx = -0.32 meV, g = 90 /ueV, ja = 176 fieV, 7^ = 133 /ieV, Pa = 0.25 7^, = 
5.3-10"'^ jx- At these quite low pump intensities the Bose-like spectrum ([9]) differs only slightly from 
that calculated including 2-photon states (dashed line). The dotted line (red) emission spectrum 
calculated according to Ref. 3] differs significantly, (b): PL spectrum (thin continuous line) from 
the data reported by Reithmaier et al. together with the corresponding fit (red). Dashed and 
dotted lines (blue and green) : PL spectra with the same parameters as the fit but with = and 
Px = respectively. 



scenario they depend upon. According to Eq. ^ it is not possible that a system in the weak 
cou,l^, enters SC tkants to „ eontrast to the results o. Re!, ft 

Although at low pump intensities, our approach and that of Ref. [7| essentially repre- 
sent models of a linear Bose-like dynamics of two coupled harmonic oscillators, nontrivial 
differences can be appreciated. Indeed, Fig. [1^ displays the PL spectrum (continuous line) 
calculated from Eq. ([9]) compared with the result obtained from the model of Ref. |3] for the 
same situation (see caption of Fig. [T]for details). At these quite low pump intensities our 
Bose-like spectrum dH]) differs only slightly from that calculated including 2-photon states 
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(dashed line) and the inclusion of additional photon states does not produce appreciable 
modifications of the spectrum. On the contrary the emission spectrum calculated according 
to Ref. differs significantly (dash-dotted line). The difference originates mainly from the 
effective reduction of the damping rates due to the incoherent pump [3], a feature absent in 
our theoretical description. The model presented in Ref. j3] reproduces the emission spectra 
measured by Reithmaier et al. ^]. Although the resulting fit is in excellent agreement with 
the experimental data, some key fitted parameters (e.g. the pumping rates) appears to be 
in contrast with the experimental results. In particular, as can be inferred from Figs. 3 and 
4 of Ref. ^ , the measured PL intensities display a very small variation as a function of tem- 
perature. On the contrary they estimate pumping rates with a variation larger than a factor 
four [l^ . In addition, Reithmeier et al. provides a quite accurate estimate of the mean cav- 
ity photon-number (see Additional Information of Ref. ^) which is one order of magnitude 
lower than the result of the above montioned fit. This disagreement strongly supports the 
absence of any dependence of the damping rates on the pumping rates in agreement with 
the results of our model. Indeed, all the above experimental features are well reproduced by 
our master equation ([T]). 

Despite of its analytical simple form, eq. iQ shows a nontrivial dependence of the system 
on the ratio between the cavity and dot pumping rates. Fig. [TJ3 displays one PL spectrum 
(thin continuous line) from the data reported by Reithmaier et al. together with the corre- 
sponding fit we obtain using Eq. ([9]) (thicker continuous line). Beside the nearly coincidence 
of the two curves in this situation, Eq. ([9]) is also in very good agreement with the experimen- 
tal spectra of Ref. jj] at different detunings. As absolute information on the PL intensities is 
commonly not available from this kind of experiments, we cannot obtain absolute values for 
the pumping rates. Instead we gather from the fit the ratio Pa/ Px = 0.86. We also obtain 
g = 76 meV, 7^ = 100 /ieV, = 35 fieV. In order to evidence the impact that the pumping 
mechanism can have on SC emission lines even at very small detuning, we plotted in Fig. [lb 
two spectra with the same parameters as the fit but with P^ = (dotted line) and Pa = 
(dashed line). These plots show that, close to the SC threshold, incoherent pumping of the 
exciton level can even hinder the line splitting. This is the case of the system investigated 
by Reithmaier et al. This result put forward the importance of adressing the correct 
mixed quantum steady state in order to describe the experimental results. 

Another situation for a similar system has been recently reported in Ref. Some first 
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FIG. 2: Color online, (a) PL spectra calculated from Eq. ([9] at different detunings reproducing the 
experimental observation of the SC of a QD-cavity system reported in Ref. (b) continuous line: 
one spectrum of panel a (the one with the lowest cavity energy) compared with the corresponding 
spectrum obtained using Pa = (dotted line) and with the corresponding absorption spectrum 
(dashed line). 



results based on the present multi-photon model ([T]) are depicted in Fig. [2^. It displays 
emission spectra at different detunings reproducing one experimental observation of the SC 
of a QD-cavity system reported in Ref. Although the parameters have been chosen 
in order to fit just one of the spectra, our model provides results which agree very well 
[isl with all the spectra presented in Fig. (2b) of Ref. changing only the cavity energy. 
We can thus unravel the ratio between the cavity-mode and the exciton pumping in this 
kind of experiments. In this case we obtain Pa/ Pb = 0.13. Fig. [2)d clearly underlines the 
importance of taking into account the correct quantum steady state resulting from the 
interplay of pumping and decay (of both the cavity mode and the exciton level) in order to 
describe properly the SC experiments. It shows one spectrum of Fig. [2|i (A = 0.1 meV) 
(continuous line) together with the emission spectrum obtained neglecting the pumping Pa 
of the cavity mode and the absorption spectrum under a coherent pump. As can be clearly 
seen without the proper Pa pumping mechanism, the calculated results differ significantly 
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from the experiment. 




FIG. 3: Transition from Rabi splitting to single-emitter laser when increasing the pumping rate. 



Our model is also well suited to describe situations well beyond that of linear or lowest 
order nonlinear effects. Fig. [3] displays emission spectra obtained at different pumping 
intensities. These results have been obtained after an exact calculation of the SC quantum 
dynamics achieved truncating the photon number-states only when including larger numbers 
(and hence dealing with a larger system of equations of motion) does not produce any change 
in the emission spectrum. We used the same broadenings of Fig. [21 a coupling g = 76 fieV 
and we set Pa = and A = 0. Increasing the excitation a reduction of the Rabi splitting 
due to saturation effects is clearly observable. The upper panel show that as soon as a 
small popoulation inversion {rix > 0.5) is reached an important narrowing of the linewidth 
appears indicating a thresholdless lasing behaviour. It seems worth underlining that, within 
our model, both the linewidth narrowing and the build up of lasing are purely fermionic 
effects not to be confused with the narrowing due to the incoherent pumping found in the 
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model of Ref. [7| (see also [5|). 

In conclusion, we have given a simple theoretical framework with a reliable predicting 
character for the analysis of photoluminescence properties under incoherent excitation of 
single quantum dot microcavity devices in a steady state maintained by a continuous in- 
coherent pumping. On the contrary to currently exploited approaches in the literature, 
we constructed a completely positive quantum master equation able to provide physically 
sensible results free from artifacts do to the reduced phenomenological description. Remark- 
ably, our theoretical framework represents a simple and flexible tool naturally suitable for 
the analysis of the impact of the interplay between the two feeding mechanisms (i.e. Pa 
and Px) onto the quantum statistic of the emitted quanta at different excitation intensities 
(under current development). At low pumping rates we obtained an analytical expression 
for the emission spectrum which directly shows the very different impact of the two feed- 
ing mechanisms on the spectra. Our model showed a great flexibility proposing simple but 
sensible descriptions of various experiments in the strong-coupling regime. It provides es- 
sentially non-perturbative results which can be exploited for the efficient modeling of future 
SC optoelectronic devices. 

We thanks Fausto Rossi for helpful suggestions and discussions on Lindblad master equa- 
tions and completely positive open dynamics. 
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